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Abstract 

Diffractive photoproduction of D* ± (2010) mesons was measured with the ZEUS 
detector at the ep collider HERA, using an integrated luminosity of 78.6 pb~\ 
The D* mesons were reconstructed in the kinematic range: transverse momen- 
tum pt(D*) > 1.9 GeV and pseudorapidity \r](D*)\ < 1.6, using the decay 
D* + — > D°tt^ followed by D° — > K~ir + (+c.c). Diffractive events were identi- 
fied by a large gap in pseudorapidity between the produced hadronic state and 
the outgoing proton. Cross sections are reported for photon-proton centre-of- 
mass energies in the range 130 < W < 300 GeV and for photon virtualities 
Q 2 < 1 GeV 2 , in two ranges of the Pomeron fractional momentum Xp < 0.035 
and xp < 0.01. The relative contribution of diffractive events to the inclusive 
D* ± (2010) photoproduction cross section is about 6%. The data are in agree- 
ment with perturbative QCD calculations based on various parameterisations of 
diffractive parton distribution functions. The results are consistent with diffrac- 
tive QCD factorisation. 
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1 Introduction 



In diffractive electron-proton scattering, the proton loses a small fraction of its 
energy and either emerges from the scattering intact, ep — ^ eXp, or dissociates into 
a low-mass state N, ep — > eXN. A large gap in rapidity separates the hadronic-state X 
with invariant-mass Mx and the final-state proton (or N). 

In the framework of Regge phenomenology [1], diffractive interactions are ascribed to 
the exchange of a trajectory with vacuum quantum numbers, the Pomeron trajectory. 
In quantum chromodynamics (QCD), the diffractive factorisation theorem [2-5] states 
that the diffractive cross section, in the presence of a hard scale, can be expressed as 
the convolution of universal partonic cross sections and a specific type of parton distri- 
bution function (PDF), the diffractive PDF (dPDF). Diffractive PDFs are interpreted as 
conditional probabilities to find a parton in the proton when the final state contains a 
fast forward proton. The dPDFs [6-9] have been determined from the HERA inclusive 
measurements of the diffractive structure function (Ff), defined in analogy with the pro- 
ton structure function (i^) [10], and can be used as input for calculations of different 
diffractive processes, for example at the Tevatron and LHC [11]. 

Diffractive collisions, producing hadronic-states X including a cc pair, are a particularly 
interesting component of diffractive ep interactions. The charm-quark mass provides a 
hard scale, ensuring the applicability of perturbative QCD even for small photon virtu- 
alities (photoproduction). At leading order (LO) of QCD, two types of photoproduction 
processes can be distinguished: direct and resolved photon processes. Charm production 
mainly proceeds via direct photon reactions, in which the exchanged photon partici- 
pates as a point-like particle, directly interacting with a gluon from the incoming proton 
(photon-gluon fusion, Fig. [I]). Thus, diffractive charm production is directly sensitive to 
the gluon content of the diffractive exchange. In the resolved photon processes, the photon 
behaves as a hadron-like source of partons, one of which interacts with a parton from the 
initial proton. Further interactions between partons from the photon and the proton may 
fill the rapidity gap, leading to a suppression of the observed cross sections in diffractive 
photoproduction. For example, an eikonal model [13] predicts a cross-section suppression 
by about a factor of three for diffractive resolved photoproduction at HERA. A similar 
mechanism was proposed to explain the rate of hard diffractive events at the Tevatron, 
which is lower than the expectations based on the dPDFs measured at HERA [14]. 

This paper presents a study of diffractive charm production, ep — > eD*X'p, with exchanged- 
photon virtuality Q 2 < 1 GeV 2 . The production of charm was tagged by identification 
of a D* ± (2010) meson in the final stat^j]. The measurement is based on a sample of 

1 From now on, the notation D* will be used for both D* + and D*~ . 
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events with a large gap in pseudorapidity between the proton and the produced hadronic 
system. Diffractive charm production was measured previously at HERA in deep in- 
elastic scattering (DIS) for photon virtualities above 1.5 GeV 2 [15-18]. Recently, the HI 
Collaboration has reported a measurement of diffractive charm photoproduction with 
Q 2 < 0.01 GeV 2 [16]. The measurement reported here is performed with about six times 
larger statistics and in a larger kinematic range than the HI results. 



2 Experimental set-up 

This measurement is based on the data taken with the ZEUS detector at the ep collider 
HERA in 1998-2000, when electrons or positrons of 27.5 GeV were collided with protons 
of 920 GeV. The sample used for this study corresponds to an integrated luminosity 
C = 78. 6± 1.8 pb _1 (13.6 pb" 1 and 65.1 pb^ 1 for the e~p and e + p samples, respectively^). 

A detailed description of the ZEUS detector can be found elsewhere [19]. Only a brief 
outline of the detector components most relevant to this analysis is given here. 

Charged particles are tracked in the central tracking detector (CTD) [20], which operates 
in a magnetic field of 1.43 T, provided by a thin super-conducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, arranged in 9 superlayers, covering the polar angle 
region_| 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks is 
v(Pt)/pt = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [21] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a(E)/E = 0.18/ \/E for electrons and 
a{E)/E = 0.35/ VE for hadrons, with E in GeV. 

In 1998-2000, the forward plug calorimeter (FPC) [22] was installed in the 20 x 20 cm 2 
beam hole of the FCAL with a small hole of radius 3.15 cm in the centre to accommodate 
the beam pipe. The FPC increased the forward calorimetric coverage by about one unit 
in pseudorapidity to r\ < 5. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > ejp. The 



2 From now on, the word "electron" will be used as a generic term for both electrons and positrons. 

3 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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bremsstrahlung photons were measured with a lead-scintillator calorimeter [23] placed in 
the HERA tunnel at Z = -107 m. 



3 Kinematics and reconstruction of variables 

Diffractive photoproduction in ep scattering (Fig. [TJ, 

e(e)+p(p)^e(e')+X(X)+p(p'), 

is described in terms of the four-momenta e, e' of the beam and scattered electrons, p,p' of 
the beam and scattered protons and X of the hadronic system. The following kinematic 
variables are defined: the photon virtuality, Q 2 = —q 2 , where q = e — e', the squared 
photon-proton centre-of-mass energy, W 2 = (p + q) 2 , and the fraction of the electron 
energy transferred to the proton in its rest frame, 

p ■ q W 2 

y = — - ~ — • 

p ■ e 2p ■ e 

The reaction can be considered to proceed through the interaction of the virtual photon 
with the diffractive exchange (Pomeron, IP). This process is described by the invariant 
mass, Mx, of the hadronic system X and the fraction of the proton momentum 

(g -P')-Q 
p-q 

carried by the diffractive exchange. 

The variables W, Mx and x P were reconstructed from the hadronic final state, using a 
combination of track and calorimeter information that optimises the resolution of the re- 
constructed kinematic variables. The selected tracks and calorimeter clusters are referred 
to as energy-flow objects (EFO) [24]. The Jacquet-Blondel formula [25] 

W m = ^2E P (E - P z ) 
was used to reconstruct W, where E p is the proton beam energy and 

i 

The invariant mass of the diffractively produced system was calculated from 

M?*f(?4-(?4-(?4 
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The sums in the above equations run over the energies Ei and momentum components 
P Xl P Yl and P Zl of all EFOs. 

The variable x P was reconstructed from 

x 

which is derived neglecting the photon virtuality (Q 2 ~ for the case of photoproduction), 
the square of the four-momentum transfer at the proton vertex (t — — (p — p') 2 ), and the 
mass of the proton. 

In addition, the inelasticity z(D*) was defined as 



p-q 

where p(D*) is the four-momentum of the D* meson. In the proton rest frame, z(D*) is the 
fraction of the photon energy carried by the D* meson. This variable was reconstructed as 

{ ] {E - P Z ) ' 

where (E — Pz)d* was calculated using the energy and momentum component Pz of the 
D* meson. 

The measured values of the variables W, z(D*), Mx and x P were corrected for energy 
losses in the inactive material of the ZEUS detector and for the loss of any particle down 
the beam pipe using Monte Carlo (MC) simulations. All variables were reconstructed 
with a resolution of better than 15% over the ranges considered. 



4 Theoretical predictions 
4.1 Monte Carlo simulations 

Monte Carlo simulations were used to calculate the acceptances, to evaluate correction 
factors for the selection inefficiencies and resolutions of the ZEUS detector and to estimate 
the background. 

The MC generator Rapgap 2.08/18 [26] was used to simulate diffractive photoproduction 
of D* mesons. The simulation was performed in the framework of the resolved-Pomeron 
model [12]. The cross section is proportional to the diffractive proton structure function, 
F^, which is parameterised by the product of the probability of the Pomeron emission 
(the so-called Pomeron flux factor) and the structure function of the Pomeron. The 
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parameterisation of the Pomeron flux factor by Streng and Berger [27] was used along with 
the Pomeron structure function obtained by the HI Collaboration (HlFit2 LO) [28]. The 
contribution of the sub- leading Regge trajectory (the Reggeon), which is only significant 
for x P > 0.01, was also included. 

The ep interactions were modelled using both direct and resolved photon processes. The 
MC resolved photon component, which amounts to about 35% of the total sample, is 
dominated by heavy-flavour excitation, in which a charm quark from the photon partici- 
pates in the hard scattering. To simulate resolved photon processes, the GRV-G-LO [29] 
set of photon parton densities was used. The simulation of charm production was per- 
formed with leading-order matrix elements. Higher-order QCD effects were approximated 
by parton showers, based on the leading logarithm (LL) DGLAP splitting functions [30]. 
Contributions from bottom production with subsequent decay into a final state with D* 
were also simulated. The bottom contribution, as predicted by the MC calculations, is not 
sizeable in any part of the kinematic range and corresponds to 2-3% of the total sample. 
The masses of the heavy quarks were set to m c = 1.5 GeV for charm and mb = 4.75 GeV 
for bottom. 

The MC generators Pythia 6.156 [31] and Herwig 6.301 [32] were used to model the 
non-diffractive photoproduction of the D* mesons. The CTEQ5L parameterisation [33] 
was used in both generators for the proton PDFs. 

The hadronisation process was simulated with the Lund string model [31] in the Rapgap 
and Pythia MCs, and according to a cluster hadronisation model [34] in Herwig. 

The generated Monte Carlo events were passed through the standard simulation of the 
ZEUS detector, based on the Geant 3.13 simulation program [35], and through the ZEUS 
trigger simulation package [36]. The simulated detector responses were then subjected to 
the same reconstruction and analysis procedures as the data. For the determination of 
the acceptance and correction factors, the generated Rapgap events were re-weighted 
in the variables Mx and z(D*), and the generated Pythia and Herwig events were 
re- weighted in the variables pt{D*) and rj(D*) to improve the description of the shapes 
of the measured distributions. 

4.2 NLO QCD calculations 

The cross sections for diffractive photoproduction of D* mesons were calculated at the 
next-to-leading order (NLO) in at s , the strong coupling constant, using the fixed-flavour- 
number scheme, in which only light flavours are active in the PDFs and the heavy quarks 
are generated by the hard interaction. The calculation was performed with the FMNR 
code in the double-differential mode [37,38]. The Weizsacker- Williams approximation [39] 
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was used to obtain the virtual photon spectrum for electroproduction with small photon 
virtualities. Diffractive PDFs were used instead of the conventional proton PDFs. The 
three sets of dPDFs used in the calculations were derived from NLO QCD DGLAP fits 
to the HERA data on diffractive deep inelastic scattering: the HI 2006 Fit A, Fit B [6] 
and the ZEUS LPS+charm Fit [7] diffractive PDFs. In the ZEUS LPS+charm fit, the 
diffractive DIS data were combined with the results on diffractive charm production in 
DIS [18] to better constrain the gluon contribution. The Reggeon contribution, which 
amounts to less than 2% for x P = 0.01 and grows up to ~ 15% at x F = 0.035, was not 
included. To account for the proton-dissociative contribution, present in the HI 2006 fits, 
the corresponding predictions were multiplied by the factor 0.81 [6]. 

The calculations were performed with a s (Mz) = 0.118 GeV [40] and m c = 1.45 GeV, the 
same values used in the QCD fits to the HERA data. The fraction of charm quarks 
hadronising as D* mesons was set to /(c — ■> D*) = 0.238 [41]. The Peterson pa- 
rameterisation [42] was used for the charm fragmentation with the Peterson parame- 
ter e = 0.035, obtained in an NLO fit [43] to ARGUS data [44]. The central NLO 
QCD predictions were obtained with the renormalisation and factorisation scales set to 
Rr = fiF = fi = \/ m c + 0.5 • \pt(c) + Pt(c)]. Here, Pt(c) and Pt(c) are the transverse 
momenta of the charm and anti-charm quarks. The uncertainties of the calculations 
were estimated by varying the renormalisation and factorisation scales simultaneously 
with the charm mass to [Ir = \ip = 0.5 • fi, m c — 1.25 GeV and to fiR = fj,p — 2 • fi, 
m c = 1.65 GeV and they were found to be of the order of Variations of the charm 

mass only resulted in a ±15% uncertainty. Uncertainties on the dPDFs were not included. 

The NLO predictions are given by the sum of point-like and hadron-like processes, the 
NLO analogues of the direct and resolved photon processes defined at LO. In all NLO 
calculations, the AFG-G-HO parameterisation [45] was taken for the photon PDFs. The 
hadron-like processes, in which the photon behaves as a source of light partons, contribute 
about 10% of the FMNR cross section. The dependence on the photon PDFs was checked 
by using the GRV-G-HO parameterisation [46] and was found to be negligible. It should 
be noted that the NLO diagrams in which the photon splits into a low-mass pair of c 
and c quarks, one of which interacts with a gluon from the proton, are considered as 
point-like photon processes in FMNR while they are effectively included in Rapgap as 
resolved-photon processes with heavy-flavour excitation. 

In the calculations of the inclusive D* photoproduction cross sections, the CTEQ5M 
parameterisation [33] with the default value of the QCD parameter (Aq^, d = 226 MeV) 
was taken for the PDFs of the proton. 
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5 Event selection and reconstruction of D mesons 

5.1 Event selection 

The events were selected online with a three-level trigger system [19,36]. At the first- 
and second-level triggers, data from CAL and CTD were used to select ep collisions and 
to reject non-ep backgrounds. At the third level, the full event information was available 
and at least one reconstructed D* candidate (see below) was required. The efficiency of 
the online D* reconstruction, relative to the efficiency of the offline reconstruction, was 
above 95%. 

Photoproduction events were selected offline by requiring that no scattered electron 
was identified in the CAL [47]. After correcting for detector effects, the most impor- 
tant of which were energy losses in the inactive material in front of the CAL 
and particle losses in the beam pipe [47,48], events were selected in the interval 
130 < W < 300 GeV (0.17 < y < 0.89). The lower limit was set by the trigger require- 
ments, while the upper limit was imposed to suppress any remaining events from deep 
inelastic scattering. Under these conditions, the photon virtuality is below 1 GeV 2 . The 
median Q 2 value was estimated from a MC simulation to be about 3 x 10~ 4 GeV 2 . 

5.2 Reconstruction and selection of D* ± (2010) mesons 

The Z}*(2010) mesons were reconstructed from the decay D* — > (D° — > Kir)ir s by means 
of the mass- difference method using charged tracks measured in the CTD. The n s particle 
from the D* decay is known as the "soft" pion because its momentum value is limited 
by the small difference between the masses of the D* and D° mesons. To ensure a good 
efficiency and a good momentum resolution, tracks were required to have pt > 0.12 GeV 
and to reach at least the third CTD superlayer. 

To reconstruct a D* candidate [49], two tracks of opposite charges were combined into a 
(Kir) pair forming a D° candidate. As kaons and pions were not identified, the mass of 
a charged kaon and a charged pion was assigned to each track in turn. 

Similarly, to form a "right-charge" track combination for a D* candidate, each (Kir) pair 
was combined with a third track (vr s ), which had the charged-pion mass assigned and 
charge opposite to that of the K meson in the (Kir) pair. To reduce the combinato- 
rial background, the tracks for the above combinations were selected with transverse 
momenta as follows: pr(K) > 0.5 GeV, Pt(^) > 0.5 GeV and Prijs) > 0.12 GeV. 
The Pt(tTs) cut was raised to 0.25 GeV for a data sub-sample, corresponding to an 
integrated luminosity of 16.9 ± 0.4 pb -1 , for which the reconstruction efficiency of 
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low momentum tracks was smaller due to the operating conditions of the CTD [50]. 
The Z}*-meson candidates were accepted provided the invariant- mass value M(Ktt) was 
consistent with the nominal M(D°) mass given by the PDG [40]. To take the mass res- 
olutions into account, the following requirements were applied, depending on pt(D*), 
the transverse momentum of the D* meson [51]: 

1.82 < M(Kn) < 1.91 GeV for Pt( d *) < 3 - 25 GeV > 

1.81 < M(Kix) < 1.92 GeV for 3.25 < p T (D*) < 5 GeV, 

1.80 < M(Ktt) < 1.93 GeV for 5 < p T {D*) < 8 GeV and 

1.79 < M(Ktt) < 1.94 GeV for Pt(D*) > 8 GeV. 



To suppress the combinatorial background further, the transverse momentum of the D* 
candidates was required to exceed 1.9 GeV and a cut on the ratio pr(D*) / E^ >10 ° > 0.1 
was applied. Here E^ 10 is the transverse energy measured in the CAL outside a cone 
of 9 = 10° around the forward direction. Monte Carlo studies showed that such a cut 
removes a significant fraction of the background whilst preserving most of the produced 
D* mesons. The measurements were restricted to the pseudorapidity range \r](D*)\ < 1.6, 
where the CTD acceptance is high. A clear signal was observed in the resulting mass 
difference AM = M(K7nr s ) — M(Ktt) distribution (not shown) at the nominal value. 

To determine the number of D* mesons in the signal range, 0.1435 < AM < 0.1475 GeV, 
the combinatorial background was modelled by "wrong-charge" track combinations 
and subtracted, after normalisation to the "right- charge" distribution in the range 
0.15 < AM < 0.17 GeV. A "wrong-charge" track combination for a (Kir) pair was 
defined as two tracks of the same charge with a soft pion (n s ) of the opposite charge. This 
subtraction yielded a signal of 12482 ± 208 inclusive D* mesons. 



5.3 Selection of diffractive events 

Diffractive events were identified by the presence of a large rapidity gap (LRG) be- 
tween the beam pipe, through which the scattered proton escaped detection, and the 
hadronic-system X [52]. The events with a LRG were selected by applying a cut on the 
pseudorapidity ?7 max of the most forward EFO with an energy greater than 400 MeV. 

Figure [2^ compares the measured r/ max distribution for all photoproduced D* events (after 
"wrong- charge" background subtraction) to a sum of the distributions from diffractive 
(Rapgap) and non-diffractive (Pythia) MC samples. The relative proportions of the two 
MC samples in the sum were chosen to give the best description of the shape of the data. 
The measured distribution shows two structures. The plateau at ?7 max < 3 is populated 
predominantly by diffractive events, while the peak around 77 max ~ 3.5 originates mainly 
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from non-diffractive events. To select diffractive events, while rejecting the majority of 
the non-diffractive events, the energy deposited in the FPC was required to be consistent 
with zero (E F p C < 1.5 GeV). Comparison between the ?7 max distributions of these data 
and MC events (Fig. [2b) confirms the considerable reduction of non-diffractive events 
in the sample. To further reduce the fraction of non-diffractive events, a cut r/ max < 3 
was applied. This cut ensures a gap of at least two units of pseudorapidity with respect 
to the edge of the forward calorimetric coverage provided by the FPC. A cut in r/ max 
correlates with the range of accessible x P values. The requirement ?7 max < 3 restricts the 
measurement to x P < 0.035. 

After the above selections, a signal of 458 ± 30 D* mesons was found in the AM distribu- 
tion (Fig. [3]) for diffractive photoproduction in the range x P < 0.035. In order to reduce 
the contributions from the Reggeon exchange and non-diffractive background, the selec- 
tion was also performed in the restricted range x F < 0.01, where 204 ± 20 D* mesons 
were observed. 

From the comparison between the measured and MC 77 max distributions (see above and 
Fig. [2^), normalisation factors were obtained for the diffractive and the non-diffractive MC 
samples. These normalisation factors were then used to evaluate the total and differential 
fractions (/ n( j) of residual non-diffractive events in the range ?7 max < 3 and to correct all 
the measured distributions for this background bin-by-bin. The total fraction / n( j = 3.3% 
was evaluated using the Pythia MC sample. Similar calculations were performed with 
the Herwig MC sample (total / n d = 15.5%) for the purpose of systematic uncertainty 
evaluation. 

The proton-dissociative events, ep — > eXN, can also satisfy the requirements ?7 max < 3 and 
-E'fpc < 1-5 GeV if the proton-dissociative system, N, has an invariant mass small enough 
to pass undetected through the forward beam-pipe. The fraction (/ p d) of background 
proton-dissociative events was measured previously to be f p d = 16±4(syst.)% [18], where 
the quoted uncertainty is due to the modelling and extraction procedure of the proton 
dissociation contribution. The proton dissociation contribution was assumed to be inde- 
pendent of all kinematic variables and was subtracted from all measured cross sections. 

6 Systematic uncertainties 

The cross section uncertainties for x P < 0.035 and x P < 0.01 were determined separately. 
The following sources of systematic uncertainty were taken into account (uncertainties for 
the range x P < 0.01 are given in brackets): 

• the CAL simulation uncertainty was determined by varying the CAL energy scale by 
±2% and the CAL energy resolution by ±20% of its value. The CAL first-level trigger 
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efficiencies were varied by their uncertainty. These variations resulted in a combined 
(±2.3)% uncertainty on the cross section; 

• the tracking-simulation uncertainties were obtained by varying all momenta by ±0.3% 
(magnetic field uncertainty) and by changing the track momentum and angular res- 
olutions by ^io% of their values. The systematic uncertainty due to the simulation 
of the track inefficiency [53] was found to be negligible. The variations yielded a 
combined cross-section uncertainty of l^J 

• the uncertainty in the subtraction of the combinatorial background was estimated 
by tightening separately by 2 MeV the lower and the upper boundary of the 
region in which the "wrong-charge" background was normalised. This contributed 
±0.2 (—0.5)% to the cross-section uncertainty; 

• the uncertainty in the FPC energy scale, evaluated by ±10% variations of the FPC 
energy in the MC, gave a systematic uncertainty of ^4 (—0.2)%; 

• the uncertainty in the selection of diffractive events was evaluated by varying the 
-Z?fpc cut by ±0.5 GeV, which yielded a cross-section uncertainty of t^'c, (^3)%, and 
the ?7 max cut by ±0.2, which yielded a cross-section uncertainty of j^jj (^0)%. The 
resulting uncertainty on the selection of diffractive events was ^ 1 (— D^o? 

• the uncertainty from the model dependence of the acceptance corrections was deter- 
mined by varying the re- weighting factors of the MC samples by ±20% of their values. 
The resulting cross-section uncertainty was ^4 (^3)% ] 

• the uncertainty from the model dependence of the non- diffractive event rejection was 
determined using Herwig instead of Pythia, yielding a cross-section variation of 
-11.9 (-6.8)%. 



The above systematic uncertainties were added in quadrature to determine the total 
systematic uncertainty. 

The overall normalisation uncertainties due to the luminosity measurement (±2.2%) and 
the D* and DO decay branching ratios (±2%) were not included in the total systematic 
uncertainty. The cross section uncertainty due to the subtraction of the proton dissocia- 
tion (±4.8%) is given separately. 
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7 Results 



7.1 Cross sections 

The differential cross section for ep — > eD*X'p in a given variable £ was calculated from 

da _ N D * ■ (1 - / nd ) ■ (1 - / pd ) 
d£ A-C-B-A^ 

where Nd* is the number of D* mesons observed in a bin of size A£. The overall 
acceptance was A = 13.9%. The combined D* — > (D° — > AV) 7r s decay branching ratio 
is B = 0.0257 ±0.0005 [40]. 

The cross sections for diffractive D*-meson photoproduction were measured in the kine- 
matic range Q 2 < 1 GeV 2 , 130 < W < 300 GeV (0.17 < y < 0.89), p T {D*) > 1.9 GeV, 
\r](D*)\ < 1.6 and x P < 0.035. No restriction in t was applied. The cross section, inte- 
grated over this range, is 

Vep^eD*x>p{x P < 0.035) = 1.49 ± 0.11(stat.)+°;^(syst.) ± 0.07(p.d.) nb. 

The last uncertainty is due to the subtraction of the proton-dissociative background (see 
Section I5T3"|) . 

The measurement was also repeated in the narrower range x P < 0.01, where the non- 
diffractive background admixture is smaller and the Reggeon contribution is expected 
to be negligible. The cross section integrated over the above kinematic region but for 
x P < 0.01 is 

Ve P ->eD*x' P (x F < 0.01) = 0.63 ± 0.07(stat.)l°;^(syst.) ± 0.03(p.d.) nb. 

For both x P ranges, the differential cross sections, measured as functions of the variables 
x P , M x , Pt(D*), r](D*), z(D*) and W, are presented in Tables dH and Figs. UK 

Figure H] compares the measured cross sections to the expectations from the 
resolved-Pomeron model calculated by means of the Rapgap MC program without 
re- weighting (Section 14. ip . To compare the shapes with the measured cross sections, the 
model prediction was normalised by a factor 0.34. Reasonable agreement between the 
shapes of the calculated and measured differential cross sections is observed. The relative 
contribution of resolved photon processes predicted by Rapgap increases towards forward 
r)(D*), small z(D*) and large Mx- 

Figures [5J1Z1 compare the measurements to the three sets of NLO predictions obtained 
from the FMNR calculations using the HI 2006 Fit A, Fit B and ZEUS LPS+charm Fit 
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dPDFs. The estimated calculation uncertainties (see Section H~2"l) are shown as the shaded 
band only for HI 2006 Fit A and are similar for other calculations. The uncertainties of 
the NLO QCD predictions are larger than the experimental ones in most bins. 

The NLO QCD calculations reproduce the x P differential cross section (Fig. [5]), in both 
shape and normalisation. A similar agreement between the calculations and the data 
is seen in Figs. [6] and H for the pt{D*), rj(D*), Mx and W differential cross sections 
in both ranges x P < 0.035 and x P < 0.01. The shapes of the differential distributions 
da/dz(D*) are not well reproduced by the NLO calculations. A better shape description 
of the z(D*) distributions is provided by Rapgap (Fig. H]). The agreement between 
the NLO QCD predictions and the data supports the validity of the QCD factorisation 
theorem in diffraction, implying the universality of diffractive PDFs. However, given the 
large experimental and theoretical uncertainties and the small hadron-like contribution 
expected by the NLO calculations, a suppression of the hadron-like component cannot be 
excluded. 



7.2 Fraction of D meson diffractive photoproduction 

The fraction of the diffractive to the inclusive (ep — > eD*Y) photoproduction cross sections 
for D* mesons was evaluated as 

TZ V (D*) = °~ e P~' eD * x 'p( Xlp < °- 035 ) 

In the kinematic region Q 2 < 1 GeV 2 , 130 < W < 300 GeV (0.17 < y < 0.89), 
Pt{D*) > 1.9 GeV and \rj(D*)\ < 1.6, diffractive production for x P < 0.035 contributes 

K V (D*) = 5.7 ± 0.5(stat.)+^(syst.) ± 0.3(p.d.)% 

to the inclusive D* photoproduction cross section. Systematic uncertainty partly cancel 
in this ratio. The residual systematic uncertainty is dominated by the measurement of the 
diffractive cross section. For the inclusive cross sections, the acceptance corrections were 
estimated with Herwig. The difference with respect to Pythia was used as a systematic 
check. 

This fraction IZ-p agrees with the values measured at HERA for diffractive DIS in similar 
kinematic ranges [15,17,18]. As shown in Fig. [8], TZt> is approximately independent of Q 2 . 

The differential dependences of the fraction TZ-d on pt(D*) , rj(D*), z(D*) and W are shown 
in Tables [3] and |4] and Fig. [9l Similar to the measurement in diffractive deep inelastic 
scattering [18], the fraction of the diffractive contribution decreases with increasing pt{D*) 
and rj(D*). The value of TZx> shows no strong dependence on either W or z(D*). 
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The NLO QCD predictions for !Zx> were obtained as the ratio of the diffractive cross 
section, calculated with the HI 2006 or ZEUS LPS+charm dPDFs, and the inclusive 
cross section, obtained with the CTEQ5M proton PDFs. The calculated ratios repro- 
duce the measured dependence of TZx> on the kinematic variables well both in shape and 
normalisation (Fig. [9]), supporting diffractive QCD factorisation. 

8 Conclusions 

Diffractive cross sections and their fraction of the total photoproduction cross 
section of D* ± (2010) mesons have been measured with the ZEUS detector at HERA 
using an integrated luminosity of 78.6 pb _1 . The D* mesons were reconstructed 
with pt > 1.9 GeV and \rj\ < 1.6. The measurements have been performed in the 
kinematic region Q 2 < 1 GeV 2 , 130 < W < 300 GeV (0.17 < y < 0.89), for the two 
ranges x P < 0.035 and x P < 0.01. 

The measured differential cross sections and the fraction of the inclusive photoproduction 
of D*^ mesons due to diffractive exchange have been compared to the predictions of 
NLO QCD calculations using available parameterisations of diffractive PDFs. The NLO 
predictions based on HI 2006 fits A and B as well as the ZEUS LPS+charm fit are 
consistent with the data. The measured fraction of D*^ meson photoproduction due to 
diffractive exchange is consistent with the measurements of D*^ meson production in 
diffractive deep inelastic scattering. Within the experimental uncertainties, this fraction 
shows no dependence on Q 2 and W . 

The results demonstrate that diffractive open-charm photoproduction is well described by 
the dPDF parameterisations extracted from diffractive DIS data, supporting the validity 
of diffractive QCD factorisation. However, given the large experimental and theoretical 
uncertainties and the small hadron-like contribution expected by the NLO calculations, 
a suppression of the hadron-like component cannot be excluded. 
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da/dxp ( nb) 


0.0 -r- 0.004 


51 ±11 t I 


0.004 0.007 


77 ±14 ± I 


0.007 -T- 0.010 


63 ±12 t e 


0.010 -r- 0.015 


47.7 ± 6.5 ± H 


0.015 + 0.020 


39.6 ± 8.7+ H 


0.020 -T- 0.025 


26.7 ± 8.5 t 


0.025 -T- 0.035 


27.0 ± 6.3 + 



Table 1: Differential cross section for diffractive photoproduction of D* mesons 
as a function of x P . The first column shows the bin ranges. The first and the second 
uncertainties are respectively statistical and systematic. The overall normalisation 
uncertainties due to the luminosity measurement (2.2%), the D* and D° branching 
ratios (2%) and the proton-dissociative contribution subtraction (4-8%) are not 
indicated. 



M x 


da/dMx (pb/GeV) 


(GeV) 


Xp < 0.010 


Xp < 0.035 


6 4- 13 


31.5 ±5.7 t H 


31.9 ±5.8 t tl 


13 4 20 


42.9 ± 7.4 1 H 


62.3 ±8.8 t H 


20 4 27 


13.4 ± 2.9 t H 


57.5 ± 7.5 t 7% 


27 4 34 


0.04 ± 0.74+°;™6 


36.2 ± 7.4 t H 


34 4 42 




12.2 ±3.6 t H 


42 4 52 




5.6 ±2.8 t H 



Table 2: Differential cross sections for diffractive photoproduction of D* mesons 
as a function of Mx for the two ranges x P < 0.035 andx P < 0.01. The first column 
shows the bin ranges. The first and second uncertainties are respectively statistical 
and systematic. The overall normalisation uncertainties due to the luminosity mea- 
surement (2.2%), the D* and D° branching ratios (2%) and the proton-dissociative 
contribution subtraction (4-8%) are not indicated. 
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p T (D*) 


da/dp T (D*) (pb/GeV) 


Kv(p T (D*)) 


(GeV) 


x P < 0.010 


x P < 0.035 


(%) 


1.9 4- 2.5 


443 ± 105 1 H 


1100 ± 194 1 [11 


6.4 ± 1.2 1 5;° 


2.5 4- 3.25 


308 ± 63 1 1 


596 ± 85 1 H 


6.1 ±0.9 1 gj 


3.25 4- 4.0 


149 ± 29 1 £ 


304 ± 42 ± || 


6.0 ±0.8 1 gj 


4.0 4- 5.0 


18.3 ±4.9 1 *J 


85.8 ± 13.1 1 £1 


3.5 ±0.5 1 ^ 


5.0 4- 6.0 


9.6 ±3.4 1 


28.6 ± 6.71 |o 


2.6 ±0.6 1 gj 


6.0 -4-10.0 


0.35 ± 0.35lg;gij 


5.09 ± 1.2 1 l-° 8 


2.0 ±0.5 1 gjj 


T)(D*) 


da/dr](D*) (pb) 






< 0.010 


x P < 0.035 


(%) 


-1.6 4- -1.2 


547 ± 98 1 ^ 


904 ± 162 1 HI 


9.5 ± 1.8 + °A 

— l.Z 


-1.2 4- -0.8 


250 ± 96 1 | 


614 ± 129 1 f 8 


5.6 ± 1.2 1 g; 7 


-0.8 4- -0.4 


287 ± 68 1 H 


775 ±124 1 go 


7.1 ±1.1 Igj 


-0.4 4- 0.0 


203 ± 71 1 24 


518 ± 100 1 £i 


5.8 ±1.11^ 


0.0 4- 0.4 


158 ± 45 1 i 8 3 


394 ± 78 1 g 


5.0 ± 1.0 1 gj 


0.4 4- 0.8 


95 ± 27 ^ ^ j 


191 ± 54 1 20 


2.8 ±0.8 1 gj 


0.8 4- 1.2 


55 ± 30 1 ^ 


220 ± 69 1 |6 


4.0 ± 1.3 1 g; 8 7 


1.2 4- 1.6 


24 ± 24 1 g;jj 


213 ± 65 1 | 


4.6 ± 1.6 1 ?;J 


z(D*) 


da/dz(D*) (pb) 


7^x>(^(-D*)) 




x*. < 0.010 


x*. < 0.035 


(%) 


0.0 4- 0.2 


1080 ± 191 + 7 4 


2726 ± 328 1 HI 


5.1 ±0.6 1 gj 


0.2 4- 0.4 


960 ±315 1 137 


2438 ± 470 1 384 


5-7 ±1.1155 


0.4 4- 0.6 


735 ± 121 1 59 


1717 ± 190 1 


6.8 ±0.8 1 gj 


0.6 4- 1.0 


157 ± 46 1 tt 


234 ± 74 1 | 


5.3 ± 1.7 ± 5:J 



Table 3: Differential cross sections for diffractive photoproduction of D* mesons 
for the two ranges x P < 0.035 and x P < 0.01 and diffractive fraction !Zx> of D* 
meson photoproduction as functions ofpT{D*), r](D*) and z(D*) . The first column 
shows the bin ranges. The first and second uncertainties are respectively statistical 
and systematic. The overall normalisation uncertainties due to the luminosity mea- 
surement (2.2%), the D* and D° branching ratios (2%) and the proton- dissociative 
contribution subtraction (4-8%) are not indicated. 
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w 


da/dW ( pb/ GeV) 


TZ V (W) 


(GeV) 


x F < 0.010 


x F < 0.035 


(%) 


130 -T- 160 


2.7 ± 1.3 i ° 5 5 


8.8 ± 1.9 1 il 


3.9 ±0.9 1 g;lj 


160 -r- 190 


4.3 ±0.9 t i 


12.1 ± 1.8 t ia 


5.6 ±0.9 1 ^ 


190 225 


4.5±1.2l jjjj 


10.6 ± 1.7 t li 


6.3 ± 1.1 1 ° - 6 7 


225 265 


3.2 ±0.7 t ° l 


5.9 ±1.0 ^ 


5.9 ±1.1 ij# 


265 -r- 300 


3.2 ±0.7 t o 


6.1 ±i.ii a | 


6.7 ± 1.2 1 



Table 4: Differential cross section for diffractive photoproduction of D* mesons 
for the two ranges x P < 0.035 and x P < 0.01 and diffractive fraction 1Zt> of D* 
meson photoproduction as a function ofW. The first column shows the bin ranges. 
The first and second uncertainties are respectively statistical and systematic. The 
overall normalisation uncertainties due to the luminosity measurement (2.2%), the 
D* and D° branching ratios (2%) and the proton-dissociative contribution subtrac- 
tion (4-8%) are not indicated. 
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Figure 1: Example of charm production in diffractive ep scattering: boson-gluon 
fusion in the resolved-Pomeron model [12]. 
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Figure 2: Comparison of the measured r\ max distribution (dots) with the sum 
(solid histograms, normalised to the data) of the weighted diffractive (^RAPGAP ) 
and non- diffractive (Pythia, shaded histograms) MC distributions for (a) all in- 
clusively photoproduced events with a reconstructed D* meson and (b) events with 
E FPC < 1.5 GeV . The D* mesons with p T (D*) > 1.9 GeV and \r)(D*)\ < 1.6 were 
selected in the kinematic region Q 2 < 1 GeV 2 and 130 < W < 300 GeV. The dis- 
tributions for the non- diffractive events as predicted by Herwig MC are indicated 
by the dotted histograms. 
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Figure 3: The distribution of the mass difference, AM = M(Kirn s ) — M(Kti), 
for the D*(2010) candidates (dots) with p T (D*) > 1.9 GeV and \r)(D*)\ < 1.6, 
reconstructed for Q 2 < 1 GeV 2 , 130 < W < 300 GeV and x P < 0.035. The shaded 
band shows the signal range, in which the combinatorial background (histogram) 
modelled by the wrong-charge combinations was subtracted. The signal contains 
458 ± 30 D* mesons. 



23 



ZEUS 



o i iii 
1 



» 10 

O 

T3 



10 



-2 



" " I " " I " " I " " I " " I " " 

• ZEUS 79 pb" 1 
Xjp < 0.035 

RAPGAPX0.34 
Resolved only 



23456789 10 




^ 4.5 
5 4 

I 3 
D 2.5 

2 

1.5 

1 

0.5 




p T (D )(GeV) 



1.5 -1 -0.5 0.5 1 1.5 

T1(D*) 



■ i i i i ■ 



0.02 

o 

^0.015 



0.2 0.4 0.6 0.8 1 

* 

z(D) 



q» 0.08 
O 

^ 0.06 
1 0.04 
0.02 




10 20 30 40 50 

M Y (GeV) 



T3 



0.01 



0.005 



150 200 250 300 
W (GeV) 




0.01 0.02 0.03 

x 



IP 



Figure 4: Differential cross sections (dots) for diffractive photoproduction of 
D* mesons with respect to px(D*), r](D*), z(D*), Mx, W and x P measured 
for x P < 0.035. The inner bars show the statistical errors; the outer bars corre- 
spond to the statistical and systematic uncertainties added in quadrature. The data 
are compared to the prediction o/Rapgap (solid histograms) using the H1FU2 LO 
diffractive parton distribution parameterisation. The theoretical prediction was nor- 
malised to the data. The dashed histograms show the predicted contribution from 
resolved photon processes. 
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Figure 5: Differential cross section (dots) for diffractive photoproduction of D* 
mesons, measured with respect to x F . The inner bars show the statistical errors; 
the outer bars correspond to the statistical and systematic uncertainties added in 
quadrature. The data are compared to the NLO QCD calculations (histograms) us- 
ing the HI 2006 Fit A (solid), Fit B (dashed), both multiplied by a factor of 0.81, 
and the ZEUS LPS+charm Fit (dotted) diffractive parton distribution parameter- 
isations. The shaded bands show the uncertainties coming from variations of the 
charm-quark mass and the factorisation and renormalisation scales. 
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Figure 6: Differential cross sections (dots) for diffractive photoproduction of D* 
mesons with respect to px{D*) , rj(D*), z(D*), Mx andW, measured for x P < 0.035. 
The inner bars show the statistical errors; the outer bars correspond to the statis- 
tical and systematic uncertainties added in quadrature. The data are compared to 
the NLO QCD calculations (histograms) using the HI 2006 Fit A (solid), Fit B 
(dashed), both multiplied by a factor of 0.81, and the ZEUS LPS+charm Fit (dot- 
ted) diffractive parton distribution parameterisations. The shaded bands show the 
uncertainties arising from variations of the charm-quark mass and the factorisation 
and renormalisation scales. 
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Figure 7: Differential cross sections (dots) for diffractive photoproduction of D* 
mesons with respect to pt{D*), rj(D*), z(D*), Mx and W, measured for x P < 0.01. 
The inner bars show the statistical errors; the outer bars correspond to the statis- 
tical and systematic uncertainties added in quadrature. The data are compared to 
the NLO QCD calculations (histograms) using the HI 2006 Fit A (solid). Fit B 
(dashed), both multiplied by a factor of 0.81, and the ZEUS LPS+charm Fit (dot- 
ted) diffractive parton distribution parameterisations. The shaded bands show the 
uncertainties arising from variations of the charm-quark mass and the factorisation 
and renormalisation scales. 
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Figure 8: Fractions IZv of D* meson diffractive production cross sections mea- 
sured at HERA in DIS [15, 17, 18] and photoproduction (this measurement). The 
inner bars show the statistical errors, and the outer bars correspond to the statistical 
and systematic uncertainties added in quadrature. 
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Figure 9: Fraction of D* meson diffractive photoproduction as a function of 
Pt{D*), r)(D*), z(D*) and W. The inner bars show the statistical errors; the outer 
bars correspond to the statistical and systematic uncertainties added in quadrature. 
The data are compared to the NLO QCD calculations (histograms) using the HI 
2006 Fit A (solid), Fit B (dashed), both multiplied by a factor of 0.81, and the 
ZEUS LPS+charm Fit (dotted) diffractive parton distribution parameterisations. 
The shaded bands show uncertainties arising from variations of the charm- quark 
mass and the factorisation and renormalisation scales. 
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